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Abstract 

We note the differences between the Kovchegov equation and the Balitsky- 
JIMWLK equations as methods of evaluating high energy hard scattering near 
the unitarity hmit. We attempt to simulate some of the correlations absent in 
the Kovchegov equation by introducing two boundaries rather than the single 
boundary which effectively approximates the unitarity limit guaranteed in the 
Kovchegov equation. We solve the problem of BFKL evolution in the presence 
of two boundaries and note that the resulting T-matrix now is the same in 
different frames, which was not the case in the single boundary case. The 
scaling behavior of the solution to the Kovchegov equation is apparently now 
lost. 
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1 Introduction 



Scattering at the unitarity limit and parton saturation [1] are dual descriptions of the 
same phenomenon [2] which lies at the heart of some of the most interesting parts of 
small-x hard scattering. The QCD dynamics giving the growth of the cross section 
with increasing energy, or equivalently, giving the growth of the parton density with 
decreasing x is BFKL [3,4] evolution. Ordinary QCD, or DGLAP [5-7], evolution 
also leads to increasing parton number densities, however, in order to get large gluon 
occupation numbers characteristic of gluon saturation (Color Glass Condensate [8]) 
it is necessary that parton densities, or cross sections, grow at a fixed hard scale 
and this requires BFKL evolution. Understanding QCD wavefunctions having high 
gluon occupation number and high energy hard scattering near the unitarity limit 
are two seemingly different aspects of the same phenomenon. The partonic language 
appears to be the more interesting description because it suggests one should be able 
to produce quite dense QCD matter. Indeed, in a central coUision of high energy 
heavy ions the system produced just after the collision is believed to be a system of 
nonequilibrium gluons at high density and high occupation numbers. This system 
then evolves into an equilibrated quark-gluon plasma of lower density and lower 
occupation numbers [9]. 

Of course genuine BFKL evolution is only accurate when gluon occupation num- 
bers are not too large, or equivalently, when S'-matrix elements are not too close to 
the unitarity limit. The Balitsky equation [10] is a generalization of the BFKL equa- 
tion which should be the leading term in a systematic treatment of the dynamics 
of high energy hard scattering at or near the unitarity limit. An equivalent for- 
malism was developed by Jalilian-Marian, lancu, McLerran, Leonidov and Kovner 
(JIMWLK) [8, 11-13] in terms of small- a; evolution of QCD wavefunctions in a light- 
cone gauge. The Balitsky and JIMWLK equations are coupled equations involving 
higher and higher correlations much as in the Schwinger-Dyson equations in ordi- 
nary perturbation theory, and as such are very difficult to deal with analytically. 
However, an interesting numerical calculation has recently appeared [14] and one 
can expect a better understanding of Balitsky and JIMWLK evolution to emerge 
from such calculations. 

Kovchegov [15] has suggested a somewhat simpler equation than the Balitsky or 
JIMWLK equations to deal with scattering at or near the unitarity limit. Although 
Kovchegov originally derived his equation in the context of a high energy scattering 
on a large nucleus, this equation can also be viewed as a mean field version of the 
Bahtsky equation in which higher correlations are neglected. For example the S'- 
matrix for the scattering of a state of two QCD dipoles on a highly evolved target, 
and at a definite impact parameter, is replaced by the product of the S'-matrices 
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of the individual dipoles. While the Kovchegov equation is not so complete it does 
have the advantage of being a precise nonlinear equation for a function. Many 
interesting limits of the Kovchegov equation have been understood by analytical 
methods. While incomplete, as is any mean field like approximation, the Kovchegov 
equation is likely the best equation one can write down in terms of functions which 
has built in correct unitarity limits for high energy scattering. In Sec. 3 we give a 
derivation of the Kovchegov equation, since its understanding plays a central role in 
this paper. 

In kinematical regions where the scattering of a QCD dipole on a target is far 
from its unitarity limit the Kovchegov equation reduces to the BFKL equation. At 
a given impact parameter the dipole target T-matrix {T = 1 — S) grows with energy 
until it approaches 1 at which point the nonlinear part of the Kovchegov equation 
limits the growth of T not to be larger than 1. If 2/(5 is the transverse size of 
the scattering dipole then the saturation momentum is defined to be the value of 
Q at which T is some preassigned number, say 1/2. This saturation momentum 
is a function of rapidity Y, Qs{Y), and gives the scale separating weak and strong 
scattering. The transition between weak and strong scattering is expected to be 
rapid as one varies the scale Q. In the weak scattering region the BFKL equation 
can be viewed, roughly, as a product of two factors: The first gives an exponential 
dependence in p = ln((5^/A^) while the second factor represents a diffusion in p with 
aV playing the role of time. In Rcfs. [16] and [17] it was suggested that one could get 
an approximate analytic solution to the Kovchegov equation simply by taking the 
diffusive part of the BFKL equation and imposing an absorptive boundary near ps — 
ln{Ql/A'^). This procedure was then used to evaluate in detail the y-dependence 
of Qs and the form of T(p,Y) when p > pg. These results confirmed the broad 
picture arrived at by less complete techniques earlier [1,18-20]. Recently Munier 
and Peschanski [21] have developed a more general and more powerful procedure for 
studying the properties of the Kovchegov equation near the saturation boundary, 
confirming the corrections of the absorptive boundary approach of Refs. [16, 17] 
along the way. 

However, the issue still remains to what extend the Kovchegov equation actually 
represents the Bahtsky and JIMWLK equations at or near the saturation region. 
One exact result of the Kovchegov equation which has been known for some time is 
the Levin-Tuchin formula [22] for the S'-matrix deep in the saturation region which 

states 

S{p,b,Y) e-^^''-P^'>' (1.1) 

where the constant c = —Cf{1 — Aq) / {Nc2x{Xo)) has a value which follows from the 
Kovchegov equation. Recently, the authors of Rcf. [23] have claimed that 5* deep 
in the saturation regime has the form given by (1.1) but with a constant at least 
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a factor of 2 smaller than the c which follows from the Kovchegov equation. The 
cause for this discrepancy is the lack of fluctuations in the Kovchegov equation. Even 
earlier, Salam [24] has argued on the basis of numerical simulations that fluctuations 
are generally large in BFKL evolution and so there is no reason to expect that they 
will be small near the saturation boundary. 

It is not hard to see that fluctuations are potentially important in evolution 
not too far from the saturation boundary. As an example consider the high-energy 
scattering of a dipole of size p/ at rapidity F on a dipole of size pi at zero rapidity. 
(The relation between p and the transverse size x±_ of a dipole is p = ln(4/A^a;^).) 
We may view the problem in terms of BFKL evolution with p > Ps{y) with the 
boundary near Ps{y) constituting an absorptive boundary for diffusion inherent in 
BFKL evolution. Further the evolution can be viewed as an evolution from (pj, 0) 
to an intermediate point (p, y) and then from (p, y) to (p/, Y). We expect the final 
answer [19] to be close to T{pi — > p/) ~ e~^^~'^°^^''f''''^^'^^\ In terms of the evolution 
through the intermediate point (p, y) we can write 



at least as far as the exponential factors are concerned. The T's in (1.2) are all 
defined by BFKL evolution with an absorbtive boundary in p at pi(y), near Ps(y). 
Fig. 1 shows in the p - Y plane the boundaries and the shaded saturation region 
limiting the BFKL evolution. Clearly T{pi ^ p) ~ e"^^"'*'"-"-^"''^''^-*-' so one would 
expect that values of p in (1.2) such that p — Ps{y) > Pf — Ps(^) would contribute 
little to the integral on the right hand side of (1.2). This should be the case since 
according to unitarity T(p — > p/) should not be larger than 1, however, such is not 
the case. The dominant contribution to the right-hand side of (1.2) comes from 
regions with p — Ps{y) > Pf — Ps{Y) where manifestly T(p pf) ^ 1. This occurs 
because the boundary for T(p — > p/) is matched to the boundary for T(pj — > p), as 
shown in Fig. 1, in order to get the same answer from both sides of eq. (1.2). There 
is unitarity violating evolution also in the region with p — Ps{y) < Pf — PsiX) when 
the evolution goes from one intermediate point to another before it ends at the final 
point. For (1.2) to be a correct formula and in order to guarantee unitarity respecting 
amplitudes for each path of evolution the saturation boundary for T(p ^ p/) should 
follow the value of p and thus it should not be the same as for T{pi ^ p). In this 
case, the BFKL evolution with a saturation boundary does not fulfill (1.2). 

In this work we introduce a second absorptive barrier at p = P2{y), as shown 
in Fig. 2, such that, crudely, T(p —>■ pf) < 1 when p < p2. This second boundary 
eliminates all unitarity violating evolution between the initial point (p,, 0) and the 
final point (p/, Y) when the evolution is viewed as proceeding in two steps; from the 
initial point (pi, 0) to an intermediate point (p, y) and then from the intermediate 
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(1.2) 
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Figure 1: BFKL evolution in the presence of one saturation boundary. 



point (p, y) to the final point (p/, Y). However, when the evolution is viewed as 
proceeding through two or more intermediate steps, (pi,0) — > {p,y) — > ip',y' > 
y)... — > (p/,y), then the second boundary p2{y) docs not eliminate all unitarity 
violating evolution between one intermediate point and another, i.e., T(p p') 
may become larger than 1. The remaining unitarity violating evolution cannot be 
ehminated without going back to the Balitsky or JIMWLK equations. 

The boundary pi is the boundary corresponding to saturation in the wavefunc- 
tion of the evolved dipole pi. This boundary is naturally viewed as an approximate 
way of imposing the unitarity limit given by the Kovchegov equation when the scat- 
tering is viewed as the scattering of an elementary dipole p/ on a highly evolved 
dipole Pi- However, we may also view the process as the scattering of an elementary 
dipole Pi on the evolved dipole pf. When pf > pi one calls this evolution a back- 
ward evolution, and it has previously been used [25] in studying shadowing effects in 
scattering on nuclei. The Kovchegov equation in this backward evolution is imposed 
approximately in terms of the boundary p2. There is in fact a complete symmetry 
here; the boundary pi is to pi in forward evolution as the boundary p2 is to p/ in 
backward evolution (see Fig. 2). 

The introduction of the second boundary p2 and the resulting symmetry which it 
brings with it has another benefit. The scattering amplitude is now boost invariant 
which was not the case for the single boundary case. Thus the calculation in the 
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Figure 2: BFKL evolution in the presence of two saturation boundaries. 

center of mass, say, now gives the same answer as arises when all the evolution is 
put into one or the other of the two dipoles. 

The paper is organized as follows: Sec. 2 gives some results from dipole-dipole 
scattering in the BFKL approximations. In Sec. 3 the Kovchegov equation is derived 
emphasizing the "mean field" aspects necessary for the derivation. The approximate 
solution is given in the absorptive boundary approximation. The frame-dependence 
of scattering is demonstrated. In Sec. 4 the second boundary, p2, is introduced and 
the BFKL equation is solved in the presence of the two boundaries in the diffusion 
approximation. The positions of the boundaries and the energy dependence of the 
saturation momentum are determined. In Sec. 5 the discussion of Sec. 4 is repeated 
for running coupling evolution. 

In Sec. 4 and Sec. 5 we have solved the BFKL equations, with boundaries pi and 
P2, in a diffusion approximation. Perhaps the main change from the single boundary 
case is that the scaling law which follows from the Kovchegov equation is now badly 
violated as explained at the end of Sec. 4.2. The effect of the two boundaries, as 
compared to the single boundary, is parametrically small but in practice very large. 
For example the anomalous dimension Aq becomes A^; where 



2(1 -Ao) (P2-Pi)' 
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and where p2 — pi ~ Inl/o;^ + p/ — PisiY) {pis is the saturation hne for internal 
evolution), parametrically. Indeed the resulting corrections are so large that one 
may well doubt the diffusion approximation we have used. It's possible that our 
procedure has overcorrected the deficiencies of the Kovchegov equation, although 
our second boundary only eliminates paths of evolution which manifestly violate 
unitarity. It would be interesting to do a numerical solution of BFKL evolution in 
the presence of two boundaries to see if the diffusion approximation is accurate. 

We feel that the main conclusion of this study is that correlations may be very 
important near the saturation boundary and in the scaling region. The best way to 
see if this is really true, and to see if our procedure of calculation has hit the essence 
of the physics, is to do accurate numerical simulations of the JIMWLK equation as 
already begun in Ref. [14] to see how different those results are from those of the 
Kovchegov equation and from our two boundary procedure of imposing unitarity. 



2 Dipole-Dipole scattering; general definitions 



In this section we give some general formulas and definitions concerning the dipole- 
dipole scattering which will be used throughout this work. 

Suppose we consider the scattering of a dipole of size a; on a dipole of size x' at 
relative rapidity Y. In a frame where one of the dipoles has rapidity y and the other 
has rapidity Y — y, the forward scattering amplitude is given by 

4^2^2^2 ^(^' ^1' ^) ^(^'' ^2, Y -y) Uddiri, r2) (2.1) 
where (Tdd{fi-if2) is the dipole-dipole cross section at the two gluon exchange level 



o-dd(n, ^2) = 27rQ; r< 1 + In 



(2.2) 



with r< being the smaller of ri, r2 and r> being the larger of ri, r2 while n{x, ri, y) 
is the number density of radiated dipoles of size r\ in the wavefunction of a parent 
dipole of size x in the rapidity interval y. The expression for n{x, ri, y) which obeys 
the dipole version of the BFKL equation reads for x > ri [26] 



n{x,ri,y) = 2 



X 



27ri 



exp 



^"^^X(%-(1-A)ln^ 



TT 



and for x < ri 



n{x,ri,y) 



■I 



dX 
2ni 



exp 



2aNr 



TT 
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X(A)|/-(1-A)ln4 



x^ 



(2.3) 



(2.4) 



or 



where 



n{x, ri,y) = ( ^ ) n{ri,x, y) (for a; < ri) , 



X(A) = V^(l) - ^^(A) - - A) 



(2.5) 



(2.6) 



with ■0(A) = r'(A)/r(A). The integration contour in eqs. (2.3) and (2.4) is parallel 
to the imaginary axis with < Re(A) < 1. 

The scattering amplitude in eq. (2.1) can be evaluated in an arbitrary frame 
since the result does not depend on the frame. Often the evalution of the amplitude 

becomes especially easy in the laboratory frame where one of the dipoles carries all 
of the rapidity (this dipole evolves) and the other one has zero rapidity (this dipole 
is an elementary dipole). In this case, for y = Y , the integration over r2 can be done 
because n{x\ r2, 0) = r2S{r2 — x') and the result is 



T{x,x',Y) = J n{x,ri,Y) add{ri,x') . 



(2.7) 



Inserting (2.2) and (2.3) in (2.7), also the integration over ri can be easily done. 



27rrl 



7T 



l + lnl^ =^aVM ^ 



X 



X 



12 



\\l - A)2 



(2.8) 



leading to 
T{x,x\Y) =TTa'^x'^ J 



2ni A2(l - A)2 



exp 



x(A)y- (1- A)ln — 



X' 



X 



X 



12 



(2.9) 



3 The Kovchegov equation and its limitations 

In this section we "derive" the Kovchegov equation [15] which is probably the best 
"simple" equation including nonlinear evolution in QCD. An analytic solution to the 
Kovchegov equation in the region close to the saturation regime has been recently 
achieved [21]; analytic results for the scattering amplitude in laboratory frame and 
for the saturation momentum as a function of rapidity are obtained. These results 
agree with the ones obtained before by solving the linear BFKL equation in the 
presence of a saturation boundary [16]. In this work, we show that the scattering 
amplitude in the vicinity of the saturation regime which results from the Kovchegov 
equation does not satisfy the completeness relation and it is frame-dependent. 
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3.1 The Kovchegov equation 



Consider the high-energy scattering of a quark-antiquark dipole on a target which 
may be another dipole, a hadron or a nucleus. For convenience we view the scat- 
tering process in a frame where the dipole is left moving and the target is right 
moving. Further we suppose that almost all of the relative rapidity of the dipole 
and the target, Y, is taken by the target so that the probabihty for the left-moving 
dipole to emit gluons before scattering off the target is small. In this frame an ele- 
mentary quark-antiquark dipole scatters on a highly evolved target. We denote the 
corresponding elastic scattering amplitude by S{xq,x-^,Y) where Xq and x^ are the 
transverse coordinates of the quark and antiquark of the dipole, respectivelly. Now 
we wish to know how S{xq, Xi,Y) changes when the rapidity Y is increased by a small 
amount dY. The increase dY can be viewed either as increasing the momentum of 
the target or as increasing the momentum of the elementary dipole. If the rapidity of 
the target is increased than its wavefunction evolves further. This evolution is given 
by the functional equation derived by Jalilian-Marian, lancu, McLerran, Leonidov 
and Kovner (JIMWLK) [11-13]. On the other hand, if the rapidity of the dipole is 
increased while that of the target is kept fixed, the probability for the dipole to emit 
a gluon increases (proportional to dY) due to the change dY. In the large Nc limit 
this quark-antiquark-gluon state can be viewed as a system of two dipoles - one of 
the dipoles consists of the initial quark and the antiquark part of the gluon while 
the other dipole is given by the quark part of the gluon and the initial antiquark. 
The probability for the production of a quark-antiquark-gluon state from the initial 
quark-antiquark dipole is [27] 

dP^^dx,dY^, (3.1) 

where is the transverse coordinate of the emitted gluon while Xq^, Xq2 and x^2 with 
the transverse sizes of the dipoles as shown in Fig. 3. The probability 
dP when multiphed with the ^-matrix for the two dipole state to elastically scatter 
on the target gives the change in the ^'-matrix, dS, for dipole-hadron scattering 

^S{x,„Y) = ^/ ^'^2^ [S^'\x^„x,„Y)-S{x^„Y)] . (3.2) 

Here S^'^\xq2: ^Hu^ ^) stands for the scattering of the two dipole state on the target 
(first graph on the left-hand side of Fig. 3) while the subtracted 5'(xoi,y) is the 
virtual contribiition necessary to normalize the wavefunction [27,28]. The later 
gives the scattering of a single dipole on the target because the gluon is not in the 
wavefunction of the dipole at the time of the scattering (last two graphs in Fig. 3) . 

Eq. (3.2) is part of a set of equations derived by Balitsky [10] and it also results 
from the JIMWLK functional equations [12, 13]. It is difficult to use eq. (3.2) because 
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of the unknown S^'^\xq2,Xi2:^)- The assumption that the scattering of the two 
dipole state on the target factorizes 

S^'^\XjQ2,Xi2,Y) = S(Xj^^2i^)^i-!^12->^) ) (3-3) 

which is a sort of a mean field approximation for the gluonic fields in the target, 
leads to the Kovchegov equation [15] 

B nN r 

—5(0^1, Y) = ^J ^'^2gfr [Si^^2, y)S{x,2, Y) - S{x,„ ¥)] . (3.4) 

This equation is probably the best "simple" equation taking into account unitarity 
corrections. In his original works [15] Kovchegov used the approximation (3.3) for 
the case that the target is a large nucleus. In this case, eq. (3.3) is a reasonable ap- 
proximation. The factorization (3.3) is however less obvious for an arbitrary target. 
It is the mean field approximation (or factorization) which makes the Kovchegov 
equation treat fluctuations incorrectly. In the region where S is small, it has been 
shown in Ref. [23] that the Kovchegov equation, while giving the form of the 5*- 
matrix correctly, gives the exponential factor twice as large as the result which 
emerges when flutuations are taken into account. Furthermore, in the next sub- 
sections, we show that the scattering amplitude which results from the Kovchegov 
equation is not consistent with the completeness relation and it is frame-dependent. 
This is because the evolution of the amplitude in rapidity is devided into two succes- 
sive steps in both cases and at the intermediate step fluctuations are not properly 
taken into account. 




Figure 3: Graphs corresponding to terms in the Kovchegov equation. 



Using N'{x^j, y) = 1 — S{Xj^j, y), another useful version of the Kovchegov equa- 
tion follows 

^M^oi,>^) = f5 / ^'^2^ mx,2,Y)+N{x,2,Y) 

-N{x^,,Y) - N{x^2:Y)N{x,2,Y)] . (3.5) 
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This equation becomes easy to use when the scattering is weak since the nonhn- 
ear term N{xQ2,y)N{xi2,y) can be dropped and the hnear equation remaining is 
the dipole version [27] of the Bahtsky-Fadin-Kuraev-Lipatov (BFKL) equation [3, 4] 
whose solution is known in the saddle-point approximation [27]. In the high-energy 
regime where unitarity corrections become important or S{x^j, Y) is small, eq. (3.4) 
is easier to use since the quadratic term 5'(a^2) ^)*S'(aii2) ^) can be neglected. The 
solution in this region is also known [22,29]. Also the transition region from the 
weak to the saturation regime (sometimes also called the "scaling region" ) has been 
explored; analytic expressions for the rapidity dependence of the saturation momen- 
tum and for the scattering amplitude have been extracted by the authors of Ref. [16] 
by solving the BFKL equation in the presence of a saturation boundary. These re- 
sults have been recently confirmed by Munier and Peschanki [21] who have directly 
solved the Kovchcgov equation in the scaling region. In the next section, we show 
these results and discuss their validity. 



3.2 Solution to the Kovchegov equation 



In the vicinity of the saturation regime, the amplitude for dipole-dipole scattering in 
laboratory frame emerging from the Kovchegov equation [21] or, equivalently, from 
the BFKL evolution in the presence of a saturation boundary [16], is^ 



T{x,x',Y) 



[QliY)x'^Y-'^A 



In 



+ A 



X exp 



Ql{Y) x'^ 

7,\n\l/QliY)x'') 
AaN,x"{X^)Y 



(3.6) 



Here all the values for the constants and the notations are adopted from Ref. [16]. 
The above amplitude results from the scattering amplitude (2.9) when evaluated in 
the vicinity of the boundary line QciY), as done in [16], where 



exp 



2aNc x(Ao) 

TT 1-Ao 



Y 



Q^{Y) x^ = 

[^X"(Ao)>^]^^ 

with Ao = 0.372 determined by 

(l-Ao)x'(Ao)+x(Ao) = 0. 



(3.7) 



(3.8) 



Up to a constant factor, QciY) equals the saturation momentum QgiY), which 
separates the saturation regime from the weak scattering region. In Fig.4a we show 

^The scattering amplitude shown here follows from eq.(4.1) of Ref. [16] in the vicinity of the 
saturation boundary, including the exponential function, which was neglected in eq.(4.2) there. 
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the line Qc{Y) in the Inx^/z^ — Y plane where z is generic; BFKL evolution takes 
place on the right-hand side of the line. All paths going into the saturation region 
on the left-hand side of the line are excluded. 



The scattering amplitude (3.6) is valid in the region 



A < In 



1 



« ^/AaN,x"{Xo)Y/7^ , 



Ql{Y) x'^ 

and A is an irrelevant small constant depending on the fixed coupling a. 



(3.9) 



In the vicinity of the saturation boundary, i.e., for x''^ close to l/Ql{Y), the 
scattering amplitude in (3.6) shows a scaling behaviour since it depends only on 
Ql{y) x''^] lines with constant Ql{y) x'"^ are lines of constant scattering amplitude. 

Expression (3.6) can be viewed as follows; one takes a dipole of size x at zero 
rapidity, evolves it to many smaller dipoles (their size is of order l/(5c(^)) as the 
rapidity increases up to F, and scatters then these dipoles on an elementary dipole 
of size x' . This is shown in Fig. 4a. Equivalcntly, one may view the BFKL evolution 
also in the opposite direction, as shown in Fig. 4b, in which case the small dipole x' 
evolves to larger dipoles (size of order l/(5c(^)) which then scatter on the elementary 
dipole X. In this case the scattering amplitude takes the form 



T{x',x,Y) 



Ao)^ 



1 



,Q2(y)a;2 



l-Ao 



X exp 



In g2(y) 2;2 + A 

i,W{Ql{Y)x') 



4aiVeX"(Ao)F 



, (3.10) 



where Q\{Y^ is defined by 



Ql{Y)xi 



/2 



1 



g^(F)x2 



(3.11) 



which is obvious from Figs. 4a and 4b. This additional way of viewing the scattering 
of two dipoles turns out to be useful in our further discussion, as we will see. 

It is also useful to show the expressions for the dipole number densities in the 

vicinity of the saturation boundary. They are calcuated in the same way as the 
scattering amplitudes. The dipole number density for x> x' reads 



n{x^ x\ Y) ■ 



\x 



[g2(y)^'2]i-Ao^ 



In 



1 



Ql{Y) X 



.12 



+A 



exp 



7rln2(l/02(y),,'2y 



^aN,x!'{Xo)Y J 
(3.12) 
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(a) 



(6) 



Figure 4: The saturation boundaries Qc and Qc separate the saturated regime 
(shaded areas) from the weak interaction region. In (a) dipole x evolves from y = 
up to y = y and scatters on the elementary dipole x' while in (b) dipole x' evolves 
from y = up to y = F and scatters on the elementary dipole x. 



and it is related to the scattering amplitude as 



n{x, x', Y) , 



(3.13) 



which is obvious from eqs. (2.3), (2.9) and (3.6). Eq. (3.13) explicitly shows that 
the elementary dipole x' scatters on the evolved dipolcs x. 

Similary, for the evolution from x' to x (with x' < x), one obtains 



n(x' X, Y) — —= — 

and has the relation 



l-Ao 



InQ^(r) x^ + A 



exp 



7i\n\Ql{Y)x') 



AaN,x"{Xo)Y 



2^2 



T{x', X, Y) = -a^x 



n{x', X, Y) , 



(3.14) 



(3.15) 



as can be seen from (2.4), (2.9) and (3.10). As compared to (3.13), eq. (3.15) shows 
that the elementary dipole x scatters on the evolved dipole x'. 
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3.3 Completeness relation 



3.3.1 Internal saturation boundaries 



In this section we show that the scattering amphtude in the laboratory frame (3.6) 
resuhing from the Kovchegov equation (or the BFKL evolution with a saturation 
boundary) is not consistent with the completeness relation. This becomes obvious 
if one tries to reproduce the result (3.6) by dividing the BFKL evolution in two 
successive steps in rapidity, say from y — to y — Y/2 and then from y — Y/2 to 
y = Y as shown in Fig. 5. To do this, we start with the scattering amplitude of two 
dipoles in the center-of-mass frame which is given by (2.1) with y — Y/2 



T{x,x',Y) = J 



^2 /y*-^ (jj^ ^2 

. 2 2 2 n{x,ri,Y/2) n{x',r2,Y/2) add{n,r2] 



(3.16) 



Upon inserting the dipole number density (2.3), where x > ri, and the dipole-dipole 
cross section (2.2) in (3.16), the integration over ri can be carried out as in (2.8), 
and the amplitude becomes 



2aN, Y x' 
--X(A)--(1-A)ln-^ 



rrr ' v\ 2 [(^1^2 fdX 1 

X n{x',r2,Y/2) . 
Using relation (2.5) for x' < r2, one can rewrite the amplitude as 



X 



(3.17) 



X n(r2,x',y'/2) 



2aN, 



Y 



TT Z 



X 



;i-A)ln-3 



X 



(3.18) 



which further can be expressed in terms of dipole number densities (cf. eq. (2.3)), 



TT 



T{x,x',Y) = ^aV^ 



A^(l - Ao) 



2^2 



n{x,r2,Y/2)n(r2,x',Y/2) . (3.19) 



It is at this point that we can identify the completeness relation. Requiring this 
expression to be the same as the result for the scattering amplitude in the laboratory 
frame (3.6), or equivalently (3.10), means that the following completeness relation, 

n(x, x', Y) = n{x, r2, Y/2) n(r2, x', Y/2) , (3.20) 

has to be satisfied for the dipole number densities n evaluated in the presence of 
saturation boundaries. We note in passing that one can use eq. (2.3) to show that 
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Figure 5: BFKL evolution in two successive steps 
the internal saturation boundaries Qc and Qc- 



y/2 ^ y in the presence of 



the naive BFKL evolution without saturation boundaries fulfiUes the completeness 
relation. 

Let us now check if the right-hand side (rhs) of (3.20) equals the left-hand side 
which is given by eq. (3.12). Inserting the dipole number density (3.12) with Y — > 
Y/2 in (3.20), the right-hand side becomes 

64 f dh-2 

TT 



[rhs of (3.20)] 



27rrf 



(3.21) 



In 



1 



Ql{Y/2)rl 



+ A 



exp 



7^\n\llQliY/2)rf^ 



X 



X 



12 



1-Ao 



In 



1 



+ A 



exp 



4aiV,x"(Ao)V/2 
Ti\n\l/Ql{Y/2)x'^ 



4aiV,x"(Ao)r/2 



Ql{Y/2)x'^ 

where Ql{Y/2) follows from (3.7) with Y replaced by Y/2, and Ql{Y/2) is deter- 
mined by matching the evolution at Y/2 (sec Fig. 5), 

QI{Y/2)tI^QI{Y/2)x' . (3.22) 

In equation (3.21) the second line gives the evolution of dipole x from ?/ = up to 
y — Y/2 and the third line gives the evolution of dipole r2 from F/2 up to F. It is 
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important to note that while the evolution of dipole x is limited by the saturation 
boundary Qci the evolution of dipole r2 is limited by its own (or internal) saturation 
boundary Qc as shown by the two shaded areas in Fig. 5 in order to respect unitarity. 



Using the matching (3.22) which gives 



[Ql{Y/2)rlY-'^ Ql{Y/2)x'^ 



l-Ao 



= [Qt{Y/2)x'x''Y-'° 
= [Ql{Y)x'^Y-'° 



and the definitions 
1 



77 = In 



2 ' 
2 



Q2(r/2)r 

equation (3.21) becomes 



C = ln 



1 



aN, 



Q^^{Y/2)x^x 



2^/2 ' 



T = 



'x"{Xo)Y , 



TT 



, (3.23) 



(3.24) 



[rhs of (3.20)] = -(^)mY)x''Y-'''^ (3.25) 

\X J -7-2 7-A 



xA [f] + A] exp 



'2t 



A [C - ?7 + A] exp 



2t 



The limits in the intgration over 77 follow from the positive definite dipole number 
densities (3.9). The integration over 77 is easily done and the result is 



[rhs of (3.20)] 



X 



2 (C + 2 A) exp 



8 fx^ 
■(C+A)^ 



mY)x''] 



/211-A0 1 

— ^(2r)exp 

-7-2 o 



■3C^ + 6CA + 4A2 



2r 



+ ((C + 2A)'-2r) exp 



5C^ + 12CA + 8A2 



4t 



^Erf 



C+2A 



4r 
(3.26) 



Inserting this expression in (3.19), one obtains for the scattering amplitude 

A2 r3C^ + 6CA + 4A2' 



X 



T{x,x',Y) = 
2(C + 2A)exp 



a^x^ 



Ag(l - Ao)^ 
■(C+A)^^ 



[Ql{Y)x''Y-'°^exp 

T2 



2t 



+ ((C + 2A)"-2r)exp 



5CH12CA+8A 



21 



4t 



Erf 



C+2A 



4r 
(3.27) 



This result is not the same as the result for the scattering amphtude in laboratory 
frame (3.6). Thus, the Kovchegov equation or, equivalently, the BFKL evolution in 
the presence of a saturation boundary does not satisfy the completeness relation. 
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3.3.2 Global saturation boundary 

In this section we show that if the evolution of dipole r2 from y/2 to y is hmited 

by the saturation boundary Qc as shown in Fig. 6, rather than Qc as in the previous 
section, then the scattering amplitude resulting from the Kovchegov equation obeys 
the completeness relation. However, in this case the evolution of dipole r2, for some 
values of r2, violates unitarity, i.e. T(r2, x', F/2) > 1. To clarify this point, let's 
consider the two regions: 

r2 < l/Qciy /2): The amplitude T(r2, x', F/2) is always larger than 1 since any 
evolution path from r2 to x' goes through the internal saturation region of dipole r2 
being on the Ihs of Qc as shown by the dashed-dotted hne in Fig. 6. 

r2 > l/Qc(^/2): Each path of evolution from x to r2 to x' respects unitarity. 
However, also in this region there are paths going from an intermediate point to 
another along which T(r2, x' , Y/2) becomes larger than 1 due to the same reason as 
above. For example along the dashed-dotted-dotted line in Fig. 6. 



x' = z 




Inx^lz^ 
t 

X = z 

Figure 6: BFKL evolution in two successive steps — > y/2 — > y in the presence of 
a global saturation boundary Qc- 

The successive evolution ^ y/2 and y/2 Y in the presence of a global 
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(fixed) saturation boundary Qc as shown in Fig. 6 is given by 
64 [ (fr2 



[rhs of (3.20)] = ^ y 



27rr2 



(3.28) 



X 



X 



X exp 



^)[Q^W2)r,^]^-^°A 
g^(r/2)a;'^]'"'° 



In 



X 



12 



In 



Qiixmri 

Ql{Y/2)x^ 



+ A 



exp 



n\Ti\l/Ql{Y/2)rl) 



AaN,x"{\o)Y/2 
In(Q^(y)a;'^e-^) 



4a7V,x"(Ao)F/2 



The last two hnes (muhiphed by 8/y^) give the dipole number density n(r2, x' , Y/2) 
which vanishes at x' = e~^/Qc or near the saturation boundary Q^- To obtain this 
result, one evaluates n(r2, x', Y) in the vicinity of Qc following [16] and requires it to 
vanish at x' — e~^/Qc{Y) (instead of x' — e~^/Qc{Y)) when saturation (eq.(4.1) of 
[16]) is taken into account. This requirement opens ways of evolution which violate 
unitarity as mentioned above. Note the difference between (3.28) and (3.21); in 
(3.21) n{r2,x,Y/2) vanishes at x' = e~'^/Qc- 



Using eq. (3.22) one gets 



In 

and 



Ql{y)ri -A 
QliY/2)x^ 



\n{Ql{Y/2)Tl) + 



2(1- A. 



In 



aN, 



''X"{K)Y 



IT 



A (3.29) 



ln(l/g^(F/2)a;'^)=ln(g^(F/2)r2^) + 



2(1 - Ac) 



In 



aN, 



TT 



HQl{Y)x 



I2\ 



(3.30) 

which allow us to evaluate eq. (3.28) in the vicinity of the saturation boundary 
x' ~ l/Qc{Y) in the logarithmic approximation. Using in addition also eq.(3.23), 
we get 

[rhs of (3.20)] = mY)x-] '''' y—^ ^ HqI{Y) x'^e"^) 

X rdHQl{Y/2)rl) \n\Ql{Y/2)rl) exp^ ^^AQliXmrP' 
Jo 

which after the integration over ln((5c(^/2) r|) gives 



aNcX"iXo)Y 



(3.31) 



[rhs of (3.20)] 



n \x' 



X 



r.12 



[Qliy)x'^]'-'^ ^ 
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In 



Ql{Y)x 



12 



+ A 



(3.32) 



This result agrees with the result for the left-hand side of the completeness relation 
given by (2.3) (the additional exponential function in (2.3) becomes one in the very 
vicinity of the saturation boundary which is considered here). Note, however, that 
in order to satisfy completeness one had to violate unitarity. 



3.4 Frame-dependence 

In this section we show that the result for the scattering amplitude obtained in 
the laboratory frame (3.6) is frame-dependent. For simplicity, wc consider here the 
scattering of two dipoles in the cm. frame (see Fig. 7) and show that the result 
disagrees with the result for the scattering amplitude in the laboratory frame (3.6). 




Inx^/z^ 



X = z 



Figure 7: Dipole-dipole scattering in the presence of saturation (shaded area) in the 
center-of-mass frame. 



The scattering amplitude of two dipoles in the cm. frame is given by (3.16), or 
after the integration over ri, by (2.9). With the dipole number densities (3.12) and 
(3.14), the scattering amplitude (2.9) becomes 



27r 
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X 



X 



8 A 



[Ql{Y/2)rlY-'^A 



In- 



QliY/2)rl 



l-Ao 



Ql(Y/2)rl 
lnQl{Y/2)rl + A 



+ A 



exp 



7rln\l/Ql{Y/2)rl) 



exp 



AaN,x"{Xo)Y/2 
7rln2(g2(F/2)ri) 



4Q;7V,x"(Ao)y/2 



(3.33) 



where Ql{Y/2) is given by (3.11). For the dipole of size x, the evolution proceeds 
from y = up to y = Y/2, while for the dipole of size x' evolution goes from 
y = Y — y = to y = Y/2 as shown in Fig. 7 where the shaded saturation regions 
of the two dipoles are shown. With the definitions in (3.23), one can easily see that 
the above amplitude leads to the result (3.27) which differs from the result for the 
scattering amplitude in the laboratory frame (3.6). 



4 BFKL evolution in the presence of saturation 



In this section we solve the BFKL equation in the presence of two saturation 
boundaries, pi and p2, as shown in Fig. 8. Consider the scattering process where 
dipole X evolves from y = Ouptoy = Y and then scatters on the elementary dipole 
x'. In this frame, the boundary pi is due to parton saturation in the wavefunction 
of the evolved dipole x. The second boundary p2 excludes all paths of evolution 
from any intermediate point {p,0 < y < Y) to the final point {x', Y) with p > 
P2 which manifestly violate unitarity as discussed for the case r2 < l/Qc{Y/2) in 
sec. 3.3.2. In other words, the second boundary p2 ensures unitarity respecting 
amplitudes for each path of evolution from {x, 0) to (x', Y) when the evolution 
proceeds in two steps, from the initial point {x, 0) to an intermediate point (p, y) 
with pi < p < p2 and then from the intermediate point (p, y) to the final point {x', 
Y). However, the second boundary does not eliminate unitarity violating amplitudes 
when the evolution proceeds through two or more intermediate stages; paths of 
evolution from an intermediate point to another may violate unitarity since they 
may enter the internal saturation regimes as explained in the case of r2 > l/QciY/2) 
in sec. 3.3.2. To eliminate also this remaining unitarity violating evolution one has to 
solve the Bahtsky [10] or JIMWLK [12, 13] equation which are difficult to deal with 
analytically. The introduction of the second boundary p2 has another advantage: it 
makes the scattering amplitude boost-invariant which was not the case for a single 
saturation boundary. Thus, the calculation in the center of mass, say, now gives the 
same result as the calculation in the laboratory frame of one of the two dipoles. 

We may view the scattering process also backwards, i.e., in a frame where the 
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evolved dipole x' scatters on the elementary dipole x. In this frame, the boundary 
P2 is due to parton saturation in the wavefunction of the evolved dipole x' and pi 
excludes all paths of evolution going to x that manifestly violate unitarity. 



4.1 Diffusion in the presence of two absorptive barriers 

In this section we solve the diffusion equation in the presence of two absorptive 
barriers and use the solution to calculate the scattering amplitude in the following 
section. 

Consider the diffusion equation 

|^(P.*) = ~^(P.*) (4.1) 

with the boundary conditions 

V'(pi,t)=V(P2,i) = (4.2) 

and the initial condition 

^l^{p, 0) = 5{p - p) . (4.3) 

This boundary- value problem can be solved by the method of separation of variables. 

Specifically, we will (a) find solutions ipn{p)t) = fn{p)gn{t) of (4.1) and (4.2) and (b) 
find the general solution ilj{p,t) satisfying also the initial condition (4.3) by taking 
a suitable linear combination of the functions ipn{p,t). 

(a) Inserting ipnip,t) = fn{p)gn{t) in (4.1), we obtain 

fn{p)9'n{t) = \fM9n{t) (4.4) 

or 

g'rXt) _ I'M .4 .X 

g^{t)/A Up) ■ ^ ■ ^ 

Next, we observe that the left-hand side of (4.5) is a function of t alone, while the 
right-hand side is a function of p alone. This implies that 

g'nit) fM 



gn{t)/A fn{p) 



= -A„ (4.6) 



for some constant A^. In addition, the boundary conditions (4.2) imply that /n(pi) = 
and fn{p2) = 0. Hence ipn{p,t) is a solution of (4.1) if 

fM + An/n(p) = , fniPl) = UP2) = (4.7) 
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and 

g':^{t) + \\n9n{t) = . (4.8) 

The boundary-value problem (4.7) has a solution only if = n^7r^/(p2 — Pi)^, and 
in this case, 

fM=siJ"-^^^) (4.9) 

V P2- Pi J 

for every positive integer n — 1,2... and p\ < p < pi- Equation (4.8), in turn, 
implies that 

g^{t) = e~(^^=^^ . (4.10) 

Hence, 

^„(p,t) = c„sin { ^'^^P-P^) \ e-(^l (4.11) 

for every constant c„ to be fixed by the initial condition (4.3). 
(b) The linear combination 

^{p,t) = Vc^sin ( ^<P-pA e-a&i (4.12) 

^ \ P2-Pl J 

satifies the boundary conditions (4.1). To satisfy the initial condition (4.3), we must 
choose the constants c„ such that 

^(p, 0) = c„sin Mp-P^A = Sip - p) . (4.13) 

on the interval pi < p, p < p2. Using Fourier analysis, the constants c„ are 
Cn — / dp d{p — p) sm — ^ — sm — ^ 

P2 - Pl Jp, \ P2- Pi J P2- Pi \ P2-P1 

(4.14) 

Thus, 

^(p, t)=y sm — ^— — — sm — ^— — ^ e ("2-^1)^ * (415) 

^P2-Pl V P2-P1 ) \ P2-P1 J 

is the general solution. 
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4.2 BFKL evolution between two saturation boundaries 



Here we wish to evaluate the scattering amphtude in the laboratory frame (2.9) in 
the region between the weak and the saturation regimes in analogy to Ref. [16]. 
Hence, we define a particular line Q^O^) in the ln{x'^/x'^) — Y plane by the two 
conditions 

'^'"\'{X,)Y + ln{Qlx')=0 (4.16) 



TT 



and 



2aNr 



71 



aNr 



X"{X,)Y . 



(4.17) 



Eq. (4.16) is the saddle-point condition while eq. (4.17) is chosen so that the scat- 
tering amplitude becomes a constant along the line Qd{Y), as we will see. p2 and 
pi denote the two boundaries hmiting the BFKL evolution as shown in Fig. 8. The 
solution to (4.16) and (4.17) is 

,2 



:i - A,)x'(Arf) + xiXd) 



TT 



2(P2 - Pi) 



and 



QliY)x'^exp 



2aN, xiXd) 



Y 1- 



X"(A. 



2(P2 - PiY xiXd) 



(4.18) 



(4.19) 



TT 1 - Ad 

Eq. (4.18) determines which is a constant for fixed Y, x and x' or fixed p2 — p\- 
For large p2 — pi, Xd approaches Aq = 0.372 which is given by (3.8) since 

^2 



Ad — Ao 



1 



(4.20) 



2(p2-pi)n-Ao • 

Now we evaluate the scattering amplitude (2.9) in a region around (5d(^) by the 
saddle-point method, when oY is large. This means that we expand the exponent 
in (2.9) around Ad (since Ad satisfies the saddle-point condition) and then use (4.16) 
and (4.17). One obtains 



T{x,x',Y) 



9 

7i(\'.r' 



{Ql{Y)x''y-'^exp 



xiii - x,y 



{p2 - Plf TT 



x"(Ad)F 



X 



dX 



2m 



-. exp 



TT 



x"(Ad)r(A-Ad)' + (A-Ad)ln 



Ql{Y)x>\ 



(4.21) 



After the simple Gaussian integral, the amplitude becomes 



T(a;, x\ Y) 



Tia^x'^ 



xiii - x,r 
1 



71 



aNr, 



X 



y^AaN,x"{Xd)Y 



exp 



{p2 - PlY TT 

7,\n\llQl{Y)x'^ 
AaN,x"{Xd)Y 



X"{X,)Y 



(4.22) 
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To simplify notations, we use 



4aN,x"{Xd)Y 



TT 



and 



ip{p-pd,t) 



1 _(PzPdf 

e * 



TTt 



(4.23) 



(4.24) 



In terms of the above the amplitude becomes 
T{x,x',Y) 



1 1 



TT 



(P2 - Plf 4 



V'(p-Pd,t). (4.25) 



Here the function ip{p — Pd, t) represents the diffusive part of the amplitude since it 
satisfies the diffusion equation 



d 19^ 
— V'(p - Pd, t) = -^q^Hp - Pd, t) . 



(4.26) 



Eq.(4.24) is not the proper solution of the diffusion equation (4.26) in the presence 
of saturation. It is at this point that we impose the unitarity. To obtain the right 
solution in the presence of saturation, one has to exclude paths which go into the 
saturation region while solving (4.26). This can be realized by requiring ijj{p — pd, t) 
to vanish at the saturation boundaries. A solution to the diffusion equation with 
two saturation boundary conditions is presented in Sec. 4.1. The general solution 
for ip^p — pd,t) in the presence of saturation is given by (4.15). However, to simplify 
calculations, we will consider the scattering ampitude at large rapidities. 



t > 



4(P2 - Pif 



TV 



(4.27) 



where the function ip{p — pd,t) in (4.15) is dominated by the lowest mode n 



^s{p, t) 



sm 



p2 - pi 



7r(p-pi(0)) 

P2- P\ 



sm 



T^{p- P\) 
Pi - Pi 



(4.28) 



Note that pi < p < P2 (cf. Fig. 8), ips{pi,'t) = ips{p2,t) = 0, and p is a parameter 
which fixes the initial condition (4.3); we choose p = at t = 0. 



With (4.28) the scattering amplitude (4.25) now becomes 

- /7r(p(0)-pi(0)) 



T{x,x',Y)- 



9 



A2(l-A,)V2-Pi 



g-(l-Ad)(p-Pci(5")) sin 



P2- pi 



sm 



7i{p-pi{Y)J 



P2- pi 

(4.29) 
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X = z 



Figure 8: BFKL evolution between two shaded saturation regimes. 

In order to fix p2 and pi, we consider the amphtudes T{x, y) and T{x' , y) with 
X > z > x' and < y < Y . T{x, z, y) follows from (4.29) with p replaced by 
p' — ln(x^/z^) and Y by 



T{x,z,y)- 



2 2 



A2(l-A,)V2-Pi 



g-(l-Ad)(p'-Pd(2/)) gjj^ 



7r (p(0)-pi(0))> | ^.^^7r(p'-pi(y)) 

P2 - Pi 



P2- Pi 

(4.30) 

T{x',z,y) can also be easily obtained; one starts with (2.9), uses ln(l/(5^(y) a;'^) 
instead of ln.{Q'^{Y) x'^) in eqs. (4.16) and (4.17), and follows otherwise the steps 
needed to compute the scattering amphtude in (4.29). Following this way, one gets 



T{x',z,y) 



2 2 
7Ta X 



A^(l-A,)V2-pi 



,(1-Ad)(p'-Pd(2/)) 



;7r(p,(y)-ln(a;Vx-^))\ ^7r(p2(?/) - p'] 



P2- pi 



P2- Pl 



(4.31) 



with 
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The amplitude T{x, z, y) becomes maximum, 



d_ 
d~p' 



-T{x,z,y) 



0, 



(4.33) 



pis 



at the saturation boundary, 



(4.34) 



We emphasize that here is the saturation hne for internal evolution once x, x' 
and Y have been fixed. It is not the saturation line for our overall scattering matrix 
T{x,x',Y). We determine pi{y) by setting the amplitude at pis 

T{p' = pis,y) = 27Tx^c, (4.35) 

with c a constant of order 1 as required by unitarity. Eqs. (4.34) and (4.35) give 



Pi(y) = 



TT 



2(P2-Pi)2 x(Ad) J" 1-A, 



In 



cAi(l-A,)3(p2-pi)^e 



7rQ;2 sin ( -Trpi (0) / (p2 - Pi ) ) 

(4.36) 

Using T{x',z,y) given by (4.31) instead of T{x,z,y), one analogously obtains 



pisiy) = p2{y) - 



(4.37) 



and 



P2{y) = In 



x 



X 



12 



1- A„ 



In 



c\l{l-\,f{p,-p,fe 



Tra^ sm 



/ Tr(p2(Y)-ln{xyx'2)) \ 

\ P2-P1 ) 



2aN, x(Ad) 

TT 1 - Aw 



1 - 



TT 



x"(Arf) 



2(P2 - Pi)' X(A,) 



{Y-y). (4.38) 



Thus, 



P2 - Pi = In 



X 



X 



12 



+ 



2aN, x{^d) 
71 1 - Ad 
2 



1 - 



1- A„ 



In 



2(P2-Pi)2 x(A.). 
cAl(l-A,)3(p2-pi)^e 



or in terms of the saturation line 

P2-pi = ln(^)-pa>^) + ^ 



In 



Tra^sin (-7rpi(0)/(p2 - Pi)) 

cAg(l-A.)3(p2-pi)V 
7rQ;2sin(-7rpi(0)/(p2 - Pi)) 



(4.39) 



(4.40) 
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with 



Pis{Y) = 



2aN^ x{Xd) 

TT 1 - An 



TT 



X"(A.) 



2(P2-Pi)2 x(A,) 



In 



2„2 



cA^(l - A,)^(p2 - Pife 
7rQ;2sm(-7rpi(0)/(p2 - Pi)) 



(4.41) 



The above quantities Pi, P2, P2 — Pi, and are determined once the dipole sizes x 
and x' and their relative rapidity Y are known. 

The distance p2 — pi becomes minimum in the region ln(x^/a;'^) > PisiY) exactly 
at ln(a;^/a;'^) = Pis{Y) or x' = l/Qis{Y), in which case (4.40) reduces to 



Ap 



1- An 



In 



cAl(l-A,)3(Ap)V 
Tra^sin (— 7rpi(0)/Ap) 



(4.42) 



where Ap denotes the minimum of p2 — pi while A^ and pi(0) follow from (4.18) and 
(4.36) with P2 — pi replaced by Ap. Comparing (4.36) with (4.42), we have 

1 



-pi(0) = Ap 



1- A, 



(4.43) 



which when inserted in (4.42) leads to 



Ap 



1 



In 



cAl(l - A,)^(A p)3e^ 



(4.44) 



1-Ad 

This equation can be solved by iteration, for small a, and the solution is 

2 13 1 

Ap(a) = — In - + — In In - + 0(const.) . (4.45) 



The scaling behavior of T, as given in (4.29), and even the value of the saturation 
momentum is rather subtle. At first glance it appears that (4.29) would give an 
exponential behavior in p, a behavior identical to the scaling behavior known to 
follow from the Kovchegov equation. However, this is not the case because Pd{y) 
has a strong p-dependence also. To investigate further the p-dependence of T define 

^ - C (4.46) 

where the p and Y dependence of ( is to be determined. Differentiating (4.16) with 
respect to p gives 

d^d -C -C 



dp 2aN,x"{Xd)Y/7r 2aN,x"{Xo)Y/7r 



(4.47) 
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while differentiating (4.20) gives 

dXd TT^ d{p2 - pi) 

dp (1 - Ao)(p2 - Pi)^ dp 

Equating tfie rigfit liand side of (4.47) and (4.48) leads to 

d{p2 - Pi) ^ C(l-Ao)(p2-pi)'^ 
dp 27raN^x"iXo)y 



(4.48) 



(4.49) 



We may also use (4.36) and (4.38) to evaluate the p-dependence of pi and p2. It is 
straightforward to see that 

dp.{Y) ^^^^r 1 dip,^p,)\ 



dp \P2 - Pi dp 

and 



dp \P2 - Pi dp 

Thus 



d{p,iY) - p,{Y)) _ d{p, - Pi) _^ ^ , ^ / 1 dip2 - pi) \ _ 
dp dp \p2 - pi dp J ' 



Dropping the small term of the right hand side of (4.52) and using (4.49), one 
finds 



d{p2-pi) (1 - Ao)(p2 - Pi)' 



dp 



27raiV,x"(Ao)F 



^ (l-Ao)(p2-pi)' 
2TiaN,x"{\o)Y 



(4.53) 



While an exact solution of (4.53) looks difficult, it is easy to see the general features 
of the solution. 

When p is not too large the solution to (4.53), neglecting the final factor [ ]~^, is 

Ap 



P2 - Pi 



1 - 



(l-Ao)(p2-pi)" 

7raArcx"(Ao)y 



(4.54) 



Eq. (4.54) is valid until p2 — pi becomes on the order of (aF)^/^ at which point the 
character of the p-dependence changes and C giving 



pi(y) ^ const' 
P2{Y) p + const" . 



(4.55) 
(4.56) 
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Also, since 

d{p-Pd{Y)) ^ ^ _ . ^ d{p2 - pi) 

dp dp ^ ' ' 

one has 

p - pd{Y) = p2 - Pi + const (4.58) 
so that T, as given by (4.29), stays very close to the value 27rx^c starting at p = 

saturation momentum for the scattering of the dipole x' on the dipole x at rapidity 
Y is 

^^^^^ ^ ^^^^^ + r3i: + ^^(l-Ao)(Ap)3 

- ^-^^^ + '^(l-Ao)(Ap)3 ^^-^^^ 

with Cs a matter of choice. While (4.59) indicates, parametrically, a small percentage 
deviation from the ps that occurs using the Kovchegov equation that deviation is 
very significant. The Cg in (4.59) is determined by the exact value of T at which one 
declares saturation to have occured. Varying the value of T by, say, 2 will vary Cg 
by a similar factor which also means that T is a function of [p — ps{Y)]/ [aY j (Ap)^) 
and does not have the usual scaling behavior of the Kovchegov equation. 

Finally, it is a little curious, and perhaps unexpected, that one now has a very 
extended region in p where T is near its unitarity limit and where our formalism 
applies. This is not the case with the Kovchegov equation. As p varies from pi (F) to 
Ps(y) the interval p^ — P\ hardly changes and so the description of T hardly changes. 
What is happening is that the evolution from |/ = to y = y, for a given p, has a 
fairly well defined scale of gluons at each intermediate value of y which contribute 
to the evolution. When p changes the scale of the intermediate gluons also changes, 
but the form of the dynamics remains exactly the same. In this picture unitarity 
limits involve a rather well defined scale of gluons at each value of y, but that scale 
depends on the initial and final dipole sizes. 



4.3 Consistency checks 

In this section we show that out result for the scattering amphtude (4.29) is boost- 
invariant. If the saturation boundaries which are fixed by x, x' and Y are used for 
any evolution step in rapidity between and Y (global saturation boundaries as in 
sec. 3.3.2), then the amplitude (4.29) obeys also the completeness relation. 
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Let us start with the completeness relation (3.20). The dipole number density on 
the left-hand side of (3.20) can be calculated analogous to the scattering amplitude 
in the previous section and in the case where x > x' it is 



n{x,x',Y) = -^(^)mY)x''] 

P2 - Pi \x'y 



P2-P1 J \ P2-P1 



Using this expression, the right-hand side of (3.20) becomes 



16 r <fr2 



[rhs of (3.20)] = — I ^ (4.61) 



X 



r^l2 



^-^^ . /7r(ln(a;Vr|)-pi(y/2)) \ . MXxi^x' j x''') - p^iY)) 
P2- P\ ) \ P2-P1 



where Ql(Y/2) is given by (4.19), with Y replaced by Y/2, and Ql(Y/2) is deter- 
mined by matching the the two evolution steps at Y/2, 

Ql{Y/2)rl = Ql{Y/2)x' . (4.62) 

This equation and ^ = ln(x^/rf) allow us to write the amplitude (4.61) as 

16 [•^^ \ r/n2/A^\ l2~\^-^d 



^ I 7r(p(0)-pi(0)) \ ^.^ MHxVx'^)-p,{Y)) 

P2-P\ j \ P2-P1 

X /^(^/^)de_2/7r(e-pi(y/2))^ 



Jpi(Y/2) 2 V P2-P1 j 



' PI (y/2) 

which after the integration over ^, 

Jpi{Y/2) 2 \ p2- Pi J 4 

agrees with (4.60). Thus, our formalism satisfies the completeness relation. 

Now let us calculate the scattering amplitude of two dipolcs in the cm. frame 
starting with (3.16). To do this, in addition to n(x, r2, Y/2) given by (4.60), we need 
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also the dipole number density n(a;', r2, F/2) with x' < r2- The later is obtained in 
the same way as the scattering amplitude (4.31) in the previous section and it is 



n(x',r2,y/2) = 



P2- Pl 



QliY/2)4 



1-Ad 



(4.64) 



X sm 



n{p2{Y) - ln{xyx''))\ . f7r{p2{Y/2) - In(xVri)) 

siri 



P2-P1 



P2 - Pl 



where Qd(Y) is determined by 



/2 



QliY) x' 



(4.65) 



with Q'^iY) given in (4.19). Inserting the dipole number densities in (2.9), the 
scattering amplitude in cm. frame becomes 



Tix,x',Y) 



7ra 



16 



(Pr2 

2A2(l-A,)2(p2-pi)2y 



(4.66) 



X 



X 



^)mY/2)rlY-'^^su. 



QliY/2) r. 



sm 



p2- Pl 
nip2iY)-\n{xyx'^)) 

P2- pl 
2 U2\ 



7r{m-Pim\ . ( 7^{Hx'/rl)-p,{Y/2)) 
sm I 

P2- pl 

7r(p2(y/2)-ln(xVri)) 

olll I 

P2- pl 



With (4.65), the variable ^ = ln(a; /r2) and the equality 



sm 



7:{p2{Y) - \n{x^/x'^)) 



P2- pl 



sm 



7r(p2(r/2)-ln(a;Vri)) 



(4.67) 



sm 



P2- pl 

7r(ln(xVx'2) -pl(F))^^ [ Ti{\n{x^ /rl) - pi{Y /2))\ 



P2- pl 



sm 



P2 -Pi 



the amplitude (4.66) reads 

T(a;, x',F) - 



na^x^ 



16 



2A2(i - \,y {p2 - p,f 

/.(P(0)-P.(0)) 

V P2 - Pl 

P2{Y/2) 



mY) x'^] 



/21 



V ■ 



(4.68) 



sm 



7r(ln(a;V,T'2) -pi(F)) 



P2-P1 



X 



■ sm 



7r(e-pi(y/2)) \ 

P2 - Pl / 



After the integration over ^ (see eq. (4.63)), the resulting scattering amplitude agrees 
with the result for the scattering amplitude obtained in the laboratory frame (4.29). 
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5 BFKL evolution with running coupling and 
boundaries 



We now add a running coupling to the discussion of Sec. 4. It is difficult to use a 
global version of BFKL evolution, covering all F- values, such as given in (4.22), so we 
shall instead use a more local approach. The boundaries pi and p2 of Fig. 8 are now 
no longer straight lines, however, the problem can still be formulated in terms of a 
similar picture with curved boundaries to be determined. BFKL evolution goes from 
a; at 1/ = to a;' at 1/ = F. Between pi and p2 the normal BFKL kernel is used, but 
we demand a solution which vanishes on the boundaries. The scattering amplitude 
thus depends on x, x' and Y, the variables fixing the region of evolution. In order 
to do the actual calculation additional variables z and y are introduced, z and y are 
the dipole size and rapidity of an intermediate stage in the evolution between x and 
x'. The BFKL evolution is written in terms of z and y. In fixed coupling evolution 
these intermediate stages of evolution were not so visible because we were able to 
identify the part of the BFKL equation that corresponded to diffusion, and then it 
was possible to combine the results of global BFKL evolution, as given in (4.25), 
with the solution to the diffusion equation with boundaries, given in (4.15). We 
note, however, that T as given in (4.29) does not obey BFKL evolution, say in the 
variables x' and Y, because a variation of these variables changes the boundaries. 
The solution (4.29) represents evolution between x ai y — and x' a,t y — Y for 
fixed X, x' and Y and it is this problem that we now turn to in case the coupling 
runs. 

With the variables x, x' and Y being fixed we consider the scattering amplitude 
T(z, ?/) going from a; at y = to ^ at y, obeying the BFKL equation, and vanishing on 
the boundaries Pi{y), P2{y)- Because running of the coupling introduces an explicit 
scale it is convenient to now use variables scaled by rather than by x^. Thus in 
this section our notation is 

Pi = ln(l/AV), = In (l/AV^) , p = In (1/AV) 

pa = \n{Ql/K'), pi = In (g?/A2) , p,^\n{QllK^) (5.1) 

where Qi and Q2 are the momenta setting the boundaries, and pi and p/ represent 
the initial and final dipole sizes. 

The BFKL equation takes the form 

|. [T/a\ = X [T/a\ (5.2) 

oy TT 

with K the usual BFKL kernel. The convolution indicated in (5.2) can be written 
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dp'K{p,p')f{p') . 



more fully as 

[K X /] (p) = 
The BFKL kernel is given as 



(5.3) 



(5.4) 



and we look for a solution to (5.2) where T{p,y) vanishes when p — pi,p2- The 
running coupling is given, as usual, by 



"^^^ = Vp ■ 



(5.5) 



The problem of finding a solution to (5.2) with K given in (5.4) and where the 
solution vanishes at pi and p2 is a well defined mathematical problem, but it appears 
too difficult to solve exactly. However, following Refs. [16] and [30] we can get an 
approximate solution to the problem by expanding x(A) in (5.4) about 



X(A) ^ x(A,) + (A - A,) x'(A,) + - (A - Xa)'x"(^d) 



(5.6) 



in which case 

K{p, p') ^ e-^'-^^^" 2 



x(A<i)+x'(A,)^ + ix"(A^)^ 



5{p - p') (5.7) 



with \d to be determined much as was done in Sec. 4.2. 
Using (5.7) in (5.2) gives 



^[T/a](p,y) = ^e 



x(A<i)+x'(Arf)^ + ^x"(A<i)^ 



,(1-Ad)p 



[T/oc] {p,y). 
(5.8) 

There are three distinct regions of parameters determining how, and to what extent, 
(5.8) can be solved: 



(i) ^»1 



m 



yi/6 

P2-Pl 

yi/6 

P2- Pi 

yi/6 



< 1 
1 



In region (i) expand the 1/p coming from the running coupling in (5.8) as 

1 _ 1 jp-Pd) 

p Pd pI 



(5.9) 
(5.10) 
(5.11) 

(5.12) 
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where, as in Sec. 4.2, pd{y) = p defines a fine of constant T. Then (5.12) leads to 
diffusion in a finear potential. Condition (i) corresponds to an exponentially small 
amplitude near the boundary p2 and the problem reduces to the problem of solving 
(5.8) with a single boundary at pi which has the solution given in Ref. [16]. In 
region (ii) one may replace 1/p by 1/ pa since the range pi < p < p2 is so small that 
the finear potential has very fittle effect. In this region we can effectively use a fixed 
coupling equation locally in y, but one must take into account that pa depends on 
y. We shall shortly carry this out in detail. It appears difficult to get an analytic 
solution in region (iii), but the form of the solution is strongly constraint by the 
known solution in region (i) and (ii). 

Now we turn to the solution of (5.8) in region (ii) where p can be replaced by 
Pd- Since we expect the p-dependence of the solution to be the same, in form, to 
the fixed coupling case it is natural to write 



T/a 



sm 



P2- Pi 



f{p,y) 



(5.13) 



In addition we suppose p2 — pi is very weakly dependent on y so that (5.8) gives, 
with pi — dpi/dy, 



cos I I / + sm ' 



P2- Pi 
2N, 



irbpa 



" g-(l-Aci)(p-pd) 



P2- Pi J \ p2- Pi J dy 



(5.14) 



e(i-A.)(p-P.) sin ( <P-P^) , f 

P2- Pi 



Eq. (4.29) suggests that the p-dependence of / should be e"*^^"'*''*^^''"'''*)^''''"^, with 
the 7 Iny term due to the fact that (5.14) is an equation for T/a not for T, in which 
case (5.14) gives 



-Pi = 



2iVe 
Tlbpa 



X'(A.) 



from the cos ^^^^^z^^ part of the equation, and 



(5.15) 



^ + (1-A.)p.= ^ 
y T^opd 



X(A.)-^X"(A.)^^^-^ 



(5.16) 



coming from the sin j part of the equation. The 7/y term in (5.16) is small 

and can be dropped. Since pd defines a line of constant amplitude pd and pi must be 
exactly equal in which case (5.15) and (5.16) becomes identical to (4.16) and (4.17) 
when the form of the running coupling, (5.5), is taken into account. 
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Eqs. (5.15) and (5.16) give (4.18) along with 



dy 



Pd 



4iVc x(Ad) 
irb 1 — Art 



1 - 



TT 



X"(A. 



2 (P2 - Pi)2 X(A,) 



or 



Pd{y) 



I^N, x(Ad) 
7r6 1 - Arf 



1 - 



X"(Ao 



2(P2-Pi)^ X(A„ 



y + Cd ■ 



(5.17) 



(5.18) 



The constant q in (5.18) follows not from solving (5.17), but reflects the fact that 
we have not systematically kept terms in (5.14) and in (5.17) of size 1/y and l/^/y, 
respectively. One such term is the j/y term in (5.16). Taking a{p) ~ [b{pd + c')]~^ 
rather than a{p) ~ would also give a contribution to c^. In contrast to 

Ref. [16], here we do not have control over such terms. As before (4.18) gives A^, 
approximately as in (4.20), when p2—pi is large. Since pd has the same y-dependence 
as pi and p2 we can write 



Pl,2 



'4A^x(M 
nb 1 — Art 



1 



X"(A.) 



2 (p2 - Pir Xi^d) 



y + ci,2 



(5.19) 



where q and ci, C2 can depend on pi and pf. From (5.13) and taking 7 = 1/2 to 
guarantee the lack of y-dependence of T when p — pd is fixed, one can write 



T(p, y) = 2nx'c ^P'~P'K {l - A,) sin ( J^k^] e-(i-^.){.-Pi(.)) . (5. 

TT \ P2-P1 J 



20) 



T, as given in (5.20) and with pi given in (5.19), satisfies the BFKL equation (5.2) 
as well as the boundary conditions. In addition 



T{pi{y) + 



l-X, 



-, y) = 2t:x^c . 



(5.21) 



Eq. (5.21) is a saturation, or unitarity, condition for evolution at intermediate values 
of the rapidity. 

Eq. (5.20) is not so strong as our fixed coupling equation (4.29) which gives 
the amplitude for evolution from a given dipole of size x in terms of boundary 
lines pi and p2 which then are completely determined in terms of the evolution. 
In the present case we are able to do the evolution at large rapidities, but we are 
unable to connect to the initial dipole, or hadron. Nevertheless, there is only one 
extra unknown constant here as compared to the fixed coupling case as we now 
demonstrate. First, we note that the relationship between p2{Y) and p/ is here 
essentially as it was in the fixed coupling case. In the fixed coupling case p2 — pf 
was determined by symmetry from the determination of pi — pi which followed from 
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unitarity. However, we equally well could have determined p2 by unitarity in which 
case the condition is 



T{p = P2{y) - z r-, p/, y) = 27rr^c . 

1 — Arf 



(5.22) 



where 



-2 = 



(5.23) 



Eqs. (5.22) and (5.23) determine the curve P2{y) as the right-most value of p from 
which BFKL evolution respects the unitarity limit in going from p to pf. Thus at 
the maximum rapidity (4.38), except for some notational changes required in the 
running coupling case, remains valid. That is 



P2{Y) = P/ + — 



In 



c\l{l-\,f{p,-p,fe 
A\p2{Y)-pf) 



(5.24) 



Eq. (5.24) determines the constant C2 in (5.19) to be 



c\i{i-x,nP2-pire - 

7:^a\p2{Y) - Pf) _ 
(5.25) 

However, we are not able to determine ci. Determining ci requires evolving into 
regions where (5.10) is no longer valid and hence would require an ability to solve 
the BFKL equation keeping the running of the coupling as given in (5.13). We do 
not know how to do this. 

One can get some insight into the regions demarcated by (5.9)-(5.11) by looking 
at Pdiv) as given by (5.18) as one reduces y to the point where ^17^- is no longer so 
small. We can write 



C2 = Pf- 



'4iVcX(Ad) 
nb 1-Xd 



i-i 



2(P2-Pi)2 x(A<i) 



Y+ 



1-A. 



In 



When ^lyr- is of order 1, pd has a dominant term along with a term of size 

y^/^. But, this is exactly what happens in the region (5.9) as seen in the solution of 
Ref. [16]. Thus when ^lyr- is of order 1 there are y^^^ corrections coming both from 
the boundaries and from the linear potential. When ^^7/6^ ^ 1 the linear potential 
effects dominate the effects due to the p2 boundary. When ^T/e^ <^ 1 the boundary 
effects dominate the linear potential effects. 



1 - - 



4 ip2 - PiY x(Ad) 



(5.26) 
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